Vortons in the SO (5) model of high temperature superconductivity 
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It has been shown that superconducting vortices with antiferromagnetic cores arise within Zhang's 
SO (5) model of high temperature supercondictivity. Similar phenomena where the symmetry is not 
restored in the core of the vortex was discussed by Witten in the case of cosmic strings. It was also 
suggested that such strings can form stable vortons, which are closed loops of such vortices. Moti- 
vated by this analogy, in following we will show that loops of such vortices in the 50(5) model of high 
T c superconductivity can exist as classically stable objects, stabilized by the presence of conserved 
charges trapped on the vortex core. These objects carry angular momentum which counteracts the 
effect of the string tension that causes the loops to shrink. The existence of such quasiparticles, 
which are called vortons, could be interesting for the physics of high temperature superconductors. 
We also speculate that the phase transition between superconducting and antiferromagnetic phases 
at zero external magnetic field when the doping parameter changes is associated with vortons. 

PACS numbers: 



I. INTRODUCTION 

In the pursuit of a unified theory of high temperature 
superconductivity and antiferromagnetism, Zhang pro- 
posed the 50(5) theory of antiferromagnetism (AF) and 
d-wave superconductivity (dSC) in the high T c cuprates 
0. The order parameter for antiferromagnetism is the 
Neel vector m which is a vector under the action of the 
50(3), the group of 3-d spatial rotations. On the su- 
perconducting side, the relevant order parameter is the 
complex superconducting order parameter tp, which de- 
scribes the gap in the electron spectrum. The effec- 
tive Lagrangian for ip is invariant under the group U(l). 
The big step that Zhang originally proposed is that the 
two symmetry groups can be combined within a larger 
symmetry group, namely 5*0(5). This means that the 
three component vector m and the complex order pa- 
rameter tp can be combined to form a "superspin" vec- 
tor ft = (tpx, mi, rri2, Tris, 1P2) which transforms under the 
group 5*0(5). The presence of doping in the cuprates 
actually breaks this symmetry down to 50(3) x U(l). 
At low doping, the AF phase is favored, corresponding 
to nonzero expectation value for \fh\, = and 

(|m|) 7^ 0). As the doping is increased eventually the 
dSC phase becomes energetically favorable with (\ipj) ^ 
and (|m|) = 0. As Zhang originally discussed in 0, the 
region of intermediate doping (near the AF-dSC phase 
boundary) should be characterized by conventional su- 
perconducting vortices, but possessing antiferromagnetic 
cores. This suggestion was verified by various groups who 
looked for numerical solutions of the classical equations 
of motion for different parameters U OJ B Furthermore, 
there has been recent experimental evidence that sug- 
gests this theoretical picture may be correct 0, IS 01 H> @ • 

Our interest in this topic arises from recent work 
[Tol ITU IT3. IT^| within a completely different context, 



namely the theory of the strong interaction, QCD, at 
high baryon density [tH IT5I ] (see for a good review 
of high density QCD and a long list of references). In 
we have shown that similar vortices with nontrivial 
core structure are present within QCD at high baryon 
density for physical values of the parameters. In this 
case, the symmetry breaking parameter responsible for 
the anisotropy is the difference between the up and down 
quark masses. The effective Lagrangian which describes 
the 50(5) theory of hig h-?}, superconductivity is very 
similar to the one used in |l2( , aside from numerical con- 
stants of course. This will prove to be useful in making 
analogies between condensed matter physics and particle 
physics throughout the course of this work. We will make 
use of the results given in throughout this paper. 

In the present paper we will show that is it possible 
to have loops of dSC vortices with AF cores that are 
classically stable objects. The source of this stability is 
the presence of conserved charges trapped on the vortex 
core, leading to nonzero angular momentum. Conser- 
vation of angular momentum prevents the vortex loops 
from shrinking and eventual disappearing. This class of 
quasiparticles which generally possess nonzero angular 
momentum and charge are called vortons. The presence 
of the AF condensate is crucial, as it is what allows the 
vortons to carry angular momentum and become classi- 
caly stable quasiparticles. 

The phenomenon where a condensate forms in the core 
of a vortex, such that the vortex can form a spinning 
loop leading to classical stability, is not by any means a 
new phenomenon. This idea was considered lon g ag o in 
the context of cosmology and cosmic strings [T7lll8l ITU 

IS3, HH, H3, H3, HE l2a IS3. Our contribution here is 

the application of the previously developed (for cosmic 
strings and high density QCD) technique to the high T c 
superconductors. 

This paper is organized as follows. In section II we will 
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review the work of 0, 0, 0] where dSC vortices with 
AF cores were originally presented. A comparison will be 
made between these vortices and other vortices with non- 
trivial core structure present in high density QCD [l2| . 
Although no new results are obtained in this section, we 
believe that making a correspondence between two very 
different fields of physics is quite a useful exercise. In 
particular, applying the top olog ical (and some analyti- 
cal) arguments developed in we reproduce the result 
HI S 01 El that there is a critical value for the coupling 
constant above which the AF core is not developed. In a 
sense it is a new explanation of the phenomenon based on 
analytical (rather than numerical) calculations. Section 
III will contain our new results where we will show that 
classically stable quasiparticles called vortons are present 
within the 50(5) theory of high temperature supercon- 
ductivity. Section IV will end with concluding remarks 
and possible experimental signatures of the quasiparti- 
cles. We also formulate a conjecture that the vortons are 
responsible for the phase transition between AF and dSC 
phases at zero external magnetic field when the doping 
parameter fi is changed. 



II. VORTICES WITH NONTRIVIAL CORE 
STRUCTURE- QCD VS. HIGH TEMPERATURE 
SUPERCONDUCTIVITY 

We will begin by briefly describing the work of 12] 
where vortices with nontrivial core structure were de- 
scribed within the context of QCD at large baryon chem- 
ical potential. This will be beneficial in order to make 
analogies between condensed matter physics and particle 
physics. We should note that this is only a review section 
and no new material will be presented here. However, 
the analogy discussed below will prove to be useful for 
the analysis which follows. We will then continue with 
a review and comparison of the vortices with nontrivial 
core structure which appear within the SO(5) theory of 
high T c superconductivity Q, Q, H H . 



A. Vortices in high density QCD 

There has been a large amount of interest within the 
particle physics community on the subject of QCD, the 
theory of the strong interaction, at large baryon den- 
sity G3, Il5| (for many more references and a nice review 
see Pij|). At zero baryon density, QCD is a theory of 
quarks and gluons which are strongly coupled, such that 
confinement takes place and the observable particles are 
colorless hadrons rather than quarks and gluons. As one 
increases the baryon chemical potential the new super- 
conducting phase when the baryon symmetry is sponta- 
neously broken occurs. To explain this phenomenon, let 
us remind that in QED, the electron- electron interac- 
tion is in general repulsive, and superconductivity is a 
very subtle effect. In non-Abelian theory, QCD, simple 



one gluon exchange is always attractive in the color 3 
channel. As is well known from conventional BCS theory 
of superconductivity, an arbitrarily small attractive inter- 
action will lead to the formation of condensate of Cooper 
pairs near the Fermi surface. This is in fact what hap- 
pens in QCD at large baryon density. The ground state 
of the high density phase of QCD is characterized by a 
diquark condensate 14, 15] analogous to the condensate 
of electron Cooper pairs present in a conventional super- 
conductor. This phase of QCD is referred to as a color 
superconducting phase. The typical chemical potential 
where this phase is thought to occur (fi ~ 500 MeV, 
A ~ 100 MeV, temp T c ~ 0.6 A, where A is supercon- 
ducting gap) cannot be realized on Earth. The interest 
in this region of the QCD phase diagram is motivated by 
the fact that such densities may be realized within the 
core of compact stars, such as neutron stars p8|. 

We will not go into specific details, only state a few 
of the main features of the color superconducting phase 
of high density QCD. For the number of quark flavors 
Nf = 3 (up, down, and strange) and the number of colors 
N c = 3, the dominant part of the diquark condensate 
takes the following form [l5j : 

kLl^r ~ e a ^e abc X2, 

~ e a0 ^e abc Y c \ (1) 

where L and R represent left and right handed quarks, 
a, (3, and 7 are the flavor indices, i and j are spinor in- 
dices, a, b, and c are color indices, and XJ and V c 7 are 
complex color-flavor matrices describing the Goldstone 
bosons. This diquark condensate breaks the original sym- 
metry group SU(3) C x SU{3) L x SU(3) R xU(l) B x U(l) A 
(color gauge symmetry, left and right flavor symmetries, 
baryonic symmetry, and axial U(1)a symmetry) down 
to the diagonal subgroup SU{3) c+ n^y This diagonal 
subgroup tells us that whenever we perform an SU (3) 
color rotation, we must simultaneously perform and left 
(right) handed flavor rotation. Since color rotations are 
now linked with flavor rotations, this phase of high den- 
sity QCD with Nf — N c = 3 is referred to as the color- 
flavor locked phase (CFL) Counting the number of 
broken generators, we see that there should be 18 Gold- 
stone bosons. Of these 18 GB, 8 of them are eaten by 
the Higgs mechanism resulting in all 8 gluons acquiring a 
mass. This leaves 10 GB, an octet related to the breaking 
of SU(3) and two singlets related to U(1)b and U(1)a- 
All of these bosons (except the one related to U(1)b) 
are actually pseudo-Goldstone bosons due to the small 
explicit violation of the symmetry. In order to describe 
the low energy degrees of freedom, namely the octet of 
Goldstone bosons, one can construct the following gauge 
invariant field: 

^ = £ XgY? = exp(z^A7/.), (2) 

c 

with the ST/ (3) generators A a normalized as Tr(A a A b ) = 
25 ab and f% ~ /i 2 /(27r 2 ) being the pion decay constant 
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which can be calculated in the large fj, limit. Prior to 
the work of it was believed that the ground 

state of the CFL phase was given by S Q = diag(l, 1, 1). 
However, it was noticed that for a physical value of the 
strange quark mass (m s 3> m u ,m d ) this may not nec- 
essarily be the case. In particular, it was argued in 
that for m s > 60 MeV along with the diquark conden- 
sate ) a new K° condensation would occur and that 
S = diag(l, 1, 1) as given above would no longer repre- 
sent the true ground state of the CFL phase, but rather, 
the vacuum expectation value of the non-diagonal ele- 
ments of ~ (S2) representing the K° GB would get a 
nonzero magnitude |30|. Therefore, S would be rotated 
in some different direction in flavor space. In the physi- 
cal case where isospin symmetry is not exact (i.e. the up 
and down quarks have different masses, m d > m u ) and 
we have overall electric charge neutrality, K° condensa- 
tion occurs. The appropriate expression for E D describing 
the K° condensed ground state can be parameterized as: 



(3) 



where describes the new Goldstone mode associated 
with K condensation along with the diquark condensate 
discussed above and 9k° = V2(\K°\)/ describes the 
strength of the kaon condensation with 30] : 



the octet of Goldstone bosons 0, |3(j : 




cos t) K o 



2 ; Meff ; 

Kb 2 M 



m 2 = am u (m d + m s ), a 



3A 2 



(4) 



where m Ul m dl m s are the masses of the up, down, and 
strange quarks respectively, // is the chemical poten- 
tial, and A is the size of the color superconducting gap 
(~ 100 MeV). In order for kaon condensation (Oro 7^ 0) 
to occur, we must have thq < n e g. This leads to the 
breaking of the hypercharge U(1)y symmetry. As dis- 
cussed in 0, |3(j , the lightest degrees of freedom in the 
CFL+A" phase are the K° and K + mesons. The es- 
sential physics of these mesons can be captured with the 
following effective Lagrangian: 



c = \d i>\ 2 - 

- x(m 2 



v 2 \d^\ 2 

2\ 2 

-yj -6m 2 &T 3 $. 



(5) 



where $ = (K + , K°) is a complex doublet describing the 
K° and K + mesons, T3 is the third Pauli matrix. The 
constants f v and v have been calculated in the leading 
perturbative approximation and are given by |32l l33| : 



f = 



21-81n2 ix 2 



IS 



2tt 2 



2 1 
V =3' 



(6) 



The remaining parameters in the effective Lagrangian (J5J 
have been obtained from a more complete description of 



8m? 



-m s (m d - m u ), 



(7) 



v 2 = Moff - m l 

A 



In the case that the parameter 8m 2 in © is zero, the La- 
grangian is invariant under the symmetry group SU(2)j x 
U(1)y — * U(l) (broken down to U(l)). From topologi- 
cal arguments we know that such a Lagrangian does not 
possess vortices since the vacuum manifold is that of 3- 
sphere and therefore does not have noncontractible loops. 
In the case that Sm 2 is relatively large, then the residual 
symmetry group is U(l) x U(l) — > U(X) and the vacuum 
manifold is that of circle, leading to the formation of 
classically stable global vortex solutions. Since 8m 2 > 
then it is K° which forms the normal global strings with 
\K°(r = 0)| = and \K°(r = 00) | = rj/V2, where the 
phase varies from to 2tt as one encircles the core of the 
vortex. From these two limiting cases, it is clear that 
there should be some intermediate region that somehow 
interpolates (as a function of 8m 2 ) between the two cases. 
At some finite magnitude of 8m 2 , an instability arises 
through the condensation of A + -field inside of the core 
of the vortex. As the magnitude of md—m u decreases, the 
size of the core becomes larger and larger with nonzero 
values of both K° and K + condensates inside the core. 
Finally, at m d = m u the core of the string (with nonzero 
condensates K° and K + ) fills the entire space, in which 
case the meaning of the string is completely lost, and we 
are left with the situation when SU(2) symmetry is ex- 
act: no stable strings are possible. As discussed in [l^ . 
the K° vortex with a K + condensate on the core can be 
approximately described by the following ansatz: 



A = ^j=/(r)e* 



(8) 
(9) 



where <fi is the azimuthal angle in cylindrical coordi- 
nates, f(r) and g(r) are solutions to the classical equa- 
tions of motion obeying the boundary condition /(0) = 
0, /(oo) = 1 and g'(0) = 0, 3(00) = 0, 3(0) = 1, and a 
gives the size of the condensate at the string core (r = 0) 
determined by the parameters of the Lagrangian. The 
width of the string when symmetry is restored in the 
core, is given by where k 2 ~ (Ks ~ m o + Sm 2 ) 

(the mass scale for K°). The width of the K + conden- 
sate when symmetry is not restored can be estimated as 
1/(3, where f3 2 ~ 28m 2 is the mass difference between 
K + and K° off of the vortex. From these estimates we 
see that as 8m 2 — > 0, the width of the vortex core in- 
creases as explained qualitatively above. In order to es- 
timate the critical point where K + condensation occurs, 
we have studied the dynamics of K + in the background 
of a K° global vortex solution. This is done by substi- 
tuting the string solution (jSJ) into the energy expression 
derived from the Lagrangian JSJ and keeping only terms 
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which are quadratic order in K + . The shift in the en- 
ergy (per unit length) in the background of a K° vortex 
is then given in dimensionless variables as follows [T^ ]: 



5E 



2 2 

rj v 



d 2 fg(r)[0 + e]g(f), 



(10) 



where 



O 



1 d d . 
f df (if 



(1 



am s (m d - m u ) 
2 



COS^o)(l-/ 2 (f)), (11) 

(12) 



The problem is reduced to the analysis of the two- 
dimensional Schrodinger equation for a particle in an 
attractive potential V(f) = —(1 — cos#^o)(l — / 2 (f)) 
with /(f) being the solution of the classical equation of 
motion for K° with the boundary conditions /(0) = 
and /(oo) = 1. Such a potential is negative every- 
where and approaches zero at infinity. As is known 
from standard course in quantum mechanics [84j |. for an 
arbitrarily weak potential well there is always a nega- 
tive energy bound state in one and two spatial dimen- 
sions; in three dimensions a negative energy bound state 
may not exist. For the two dimensional case (the rele- 
vant problem in our case) the lowest energy level of the 
bound state is always negative and exponentially small 
for small A'. One should note that our specific potential 
V(f) = -(l — cos9 K o)(l — f 2 (f)) which enters (fT77f> is not 
literally the potential well, however one can always con- 
struct the potential well V such that its absolute value 
is smaller than |V(f)| everywhere, i.e. \V'\ < \V(f)\ for 
all f . For the potential well V we know that the neg- 
ative energy bound state always exists; when V is re- 
placed by V it makes the energy eigenvalue even lower. 
Therefore, the operator Ijllfl has always a negative mode 
irrespective of the local properties of function f(r). As 
a consequence, if e = then the string JSJ is an unstable 
solution of the classical equation of motion, the result we 
expected from the beginning from topological arguments. 
The instability manifests itself in the form of a nega- 
tive energy bound state solution of the corresponding 
two-dimensional Schrodinger equation i|10|) irrespective 
of the magnitudes of the parameters. The problem of de- 
termining when i^ + -condensation occurs is now reduced 
to solving the 2-d Schrodinger type equation Og = Eg. 
From the previous discussions we know that E for the 
ground state is always negative. However, to insure the 
instability with respect to i^-condensation one should 
require a relatively large negative value i.e. E + e < 0. 
It can not happen for arbitrary weak coupling constant 
~ (1 — cosOx ) when 9 K o is small. However, it does hap- 
pen for relatively large 9k°- To calculate the minimal 
critical value cr ;t when K + -condensation develops, one 
should calculate the eigenvalue £ as a function of pa- 
rameter 9 K o and solve the equation E(9 cri t) + e = 0. For 
very small coupling constant A' = (1 — cos9 K o) — > 
the bound state energy is negative and exponentially 



small, E ~ — e~~ . However, for realistic parameters of 
p. A, m s , m u , rrid the parameter e is not very small and 
we expect that in the region relevant for us the bound 
state energy E is the same order of magnitude as the po- 
tential energy ~ A'. In this case we estimate # cr it from the 
following conditions — -E(# C rit) ~ A' ~ (1 — cos0 cr ;t) ~ e 
with the result which can be parametrically represented 
as 



sin ■ 



const - 



(13) 



where we have neglected all numerical factors in order 
to explicitly demonstrate the dependence of # cr it on the 
external parameters. The limit of exact isospin symme- 
try, which corresponds to rrid — > m u when the string 
becomes unstable, can be easily understood from the ex- 
pression H13fl . Indeed, in the case that the critical pa- 
rameter # cr it — > becomes an arbitrarily small number 
the K + instability would develop for arbitrarily small 
9 K o > 0. 

The region occupied by the K + condensate at this 
point is determined by the behavior of lowest energy 
mode g at large distances, g[f — > oo) ~ exp(— Ef) such 
that a typical f ~ (m<j — m u )~ l — ► oo as expected. 



B. Vortices in the SO (5) theory of high 
temperature superconductivity 

We will now review the work of Q, i a i where 

it was shown that vortices with nontrivial core struc- 
ture similar to the ones discussed above for high den- 
sity QCD are present within the SO(5) theory of high 
T c superconductivity. The effective Lagrangian which de- 
scribes Neel vector m and dSC order parameter ip in the 
resence of zero external electromagnetic field is given by 



C = 



|(|a t m| 2 + |^| 2 )~^(|Vm| 2 + |Vm| 2 ) 



+ (g - a)\rh\ 2 - a\i>\ 2 



-b\m\ 



-&M 4 -6H 2 H 2 , 



(14) 



where we neglected the electromagnetic contribution to 
the vortex structure. Actually, one can show that the 
electromagnetic field does not change the qualitative ef- 
fects which are the subject of the present paper, and 
therefore, will be ignored in what follows. Here x is th e 
susceptibility and p = h 2 /m* is the stiffness parame- 
ter. In reality, we know that the properties of x an d 
p are different in different directions. However, when 
we discuss the topological properties of the configura- 
tion this difference can change the quantitative results 
but cannot change the qualitative picture. The Neel vec- 
tor has three spatial components m = (rni, 7712,7713) and 
the superconducting order parameter is a complex field 
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if> = tpi + iip2- The parameters of the above effective 
Lagrangian are given by: 



a < 0, b > 0, 
g = 4 X (Mc-Ai 2 )- 



P 
2\a\ 



(15) 
(16) 

(17) 



where fi is the chemical potential (or doping, not to be 
confused with the chemical potential for QCD in the pre- 
vious section), /i c is the critical chemical potential which 
defines the AF-dSC phase boundary, and £ is the coher- 
ence length. The anisotropy g is included which explic- 
itly breaks the SO (5) symmetry in the following fash- 
ion, 50(5) -> 50(3) x 17(1). If g = then the 50(5) 
symmetry is restored and the order parameters to and 
ip can be organized into a superspin order parameter 
ft = (ipi,Tni,m2,m3,i()2) which transforms in the vector 
representation of 50(5), as Zhang originally proposed in 
[lj. In the following we will consider \i > [i c and g < 
so that we are in the dSC phase and (\ip\) — \J\a\/b and 
(|m|) = in the bulk. One immediately notices that 
the form of the Lagrangian for the 50(5) theory is very 
similar to the Lagrangian used to describe K strings in 
high density QCD (JSJ in the previous section. In partic- 
ular, the key element in construction of the vortices with 
non-zero condensate in the core, the asymmetry param- 
eter, is determined by the magnitude of Sm 2 in eq. J5J. 
For the 50(5) theory it is replaced by the parameter of 
anisotropy g in eq. (|16|l . 

A nonzero vacuum expectation value for ip signals the 
onset of superconductivity and the breaking of the U(l) 
symmetry. It is well known that stable vortices can 
form since the topology of the vacuum manifold is that 
of a circle. Analogous to the situation for high den- 
sity QCD, for a certain range of the anisotropy param- 
eter g these vortices should have an antiferromagnetic 
core ((|m(r = 0)|) ^ 0). This was initially pointed 
out by Zhang pj when he introduced the 50(5) model 
and further studied in 0. 0. [| Similar to the K vor- 
tex/condensate solution given by (|8I9|) . these vortices are 
described by the following static field configurations: 




o i0 



g(r) m, 



(18) 



(19) 



where (p is the azimuthal angle in cylindrical coordinates, 
a is the parameter obeying the relation < a < 1, and 
to is an arbitrary unit vector. As before, f(r) and g(r) 
are solutions to the classical equations of motion satis- 
fying the boundary conditions /(0) = 0, /(oo) = 1 and 
g '(0) = 0, g(oo) = 0, s(0) = 1. The width of the vor- 
tex determined by the profile function / is approximately 
given by the coherence length i|17|). 6$ w £. The width of 
the condensate in the core (if it forms) is estimated to be 
of order of <L 



very large at the phase boundary. 



/j, 2 .) and becomes 



Using what we have already learned from QCD and 
the results from earlier work on these vortices [3, 0, 3 , 
we can immediately summarize the main features of these 
objects. Numerical calculations in Q confirm that as the 
anisotropy parameter \g\ is decreased, the size and width 
of the condensate in the core increases. We will support 
these numerical calculations using some analytical argu- 
ments given below. 

The free energy (per unit length) obtained from the 
Lagrangian l|14|) is: 



T 



d 2 r 



|(|5 t m| 2 + |^| 2 ) +|(|Vto| 5 



|VV| 2 )-(|a|+.9)|TO| 2 -|a||V| 5 



-&H 4 +-b\ip\ A + b\fh\ 2 \ipf 



(20) 



If anisotropy 3 = we know from topological arguments 
that this theory does not possess any vortices since the 
vacuum manifold is that of a 4-sphere and therefore does 
not have noncontractible loops. If \g\ is relatively large, 
then the residual symmetry group has a subgroup U(l) 
and the vacuum manifold is that of circle, leading to the 
formation of classically stable global vortex solution de- 
scribed in terms of ip field with a typical profile function 
when \ip(r = 0)| = and \ip(r — oo)| = \J\a\/b. From 
these two limiting cases, as discussed in the previous sub- 
section IIA, it is clear that there should be some interme- 
diate region that somehow interpolates (as a function of 
g) between the two cases. The way how this interpolation 
works is as follows (see the previous subsection where 
the physical picture is quite analogous to the present 
case). At some finite magnitude of g, an instability arises 
through the condensation of m-field inside of the core of 
the vortex. As the magnitude of g decreases, the size of 
the core becomes larger and larger with nonzero values 
of both to, and tp condensates inside the core. Finally, at 
3 = the core of the string (with nonzero condensates to 
and ip) fills the entire space, in which case the meaning 
of the string is completely lost, and we are left with the 
situation when the 50(5) symmetry is exact: no stable 
strings are possible. 

In order to estimate that critical value of the parameter 
3 = 3 cr it where an AF core forms inside the vortex, the 
same method can be applied as described for the QCD 
color superconductor in the previous section. We will use 
the following change of variables in order to express the 
free energy in terms of dimensionless variables only: 





r = £r" 



(21) 



Expanding the expression for the change in the free en- 
ergy 1)20(1 in the background of a ip vortex solution given 
by 1)18(1 and keeping only quadratic terms in to, we have: 



6JF 
T 



6 



where 



H 



1 d 



7(^x7) -(WV)), 



r' dr' dr 



(23) 
(24) 



Since we are working in the dSC phase g < the per- 
turbation e > 0. We have now cast the change in the 
free energy in the exact same form as we did for the 
QCD vortices in the previous section. The problem 
is now reduced to the analysis of the two-dimensional 
Schrodingcr equation for a particle in an attractive po- 
tential V(r') = — (1 — f 2 (r')). As before, this potential is 
negative everywhere and approaches zero at infinity. This 
means that the ground state eigenfunction Hgo = Ego 
has a negative eigenvalue E. The instability with respect 
to formation of the AF condensate in the core occurs not 
for arbitrary small negative eigenvalue E, but when the 
absolute value of |J5| is large enough to overcome the pos- 
itive contribution due to e. Therefore, we immediately 
see that an AF core forms if E + e < 0. If \g\ is greater 
than some critical value g CT n then it is not energetically 
favorable for an antiferromagnetic core to form and dSC 
vortices will possess a normal core where the symmetry 
is restored. Following the same procedure as in the QCD 
case, we have: 



\<)crit\ A X(P 2 ~ Mc) 



0.2. 



(25) 



where for numerical estimates we used the variational 
approach developed in |l9|. 

Above we have reviewed the basic properties of su- 
perconducting vortices with an antiferromagnetic core 
within the 50(5) theory of superconductivity. We should 
emphasize once more that all results presented above are 
not new and have been discussed previously from a dif- 
ferent perspective. Let us repeat the main results of this 
section once again: If g = then the dSC vortices are un- 
stable. If < \g\ < <7 cr it then an AF core will form inside 
the dSC vortices. The width of the AF core in this case 
becomes larger and larger when we approach the phase 
transition line, i.e. g — > 0. Finally, if \g\ > g CT i t then the 
dSC vortices will have a normal core when symmetry is 
restored and \m\(r = 0) = \tp\(r = 0) = 0. In what fol- 
lows we will always be working in the region of the phase 
diagram where < \g\ < g CI it and dSC vortices have an 
AF core (which will be referred to as dSC/AF vortices). 
Now we will proceed to the next section and introduce 
vortons, loops of dSC/AF vortices which are stabilized 
by angular momentum. 



III. VORTONS IN HIGH TEMPERATURE 
SUPERCONDUCTIVITY 

We will now consider the interesting possibility that 
loops of the dSC/AF vortices can exist as classically 



stable objects ( and at least metastable quantum me- 
chanically). This stability arises through a mechanism 
where topological and Noether charges can be trapped 
on the core of the vortex. Such objects, called vortons, 
have been studied extensively in the context of cosmol- 
ogy where cosmic st ring s have nontrivial core structure. 
El El 13 IM EJ El El H3. Such vortons are also 
present within high density QCD where vortices with a 
condensate trapped on the core are realized ■ 

As Davis and Shellard originally pointed out in [2(| 
I22I |23| , if one has a theory which contains vortices with 
a condensate trapped on the core then loops of such vor- 
tices can form which are stabilized by angular momentum 
alone. We will consider a large loop of string of radius 
R^> 5, where 6 is the vortex thickness, so that curvature 
effects can be neglected. The z-axis is defined along the 
length of the string, varying from to L = 2nR as one 
goes around the loop. Although we are considering a cir- 
cular loop for simplicity at the moment, we realize that 
this is probably not the relevant physical case. The re- 
sults we will discuss in this section should not depend on 
the geometry of the loop, the important point is the pres- 
ence of conserved charges which are trapped on the vor- 
tex leading to stability. In reality, the final stable config- 
uration of these vortex loops is probably a more compli- 
cated shape because of the quasi two-dimensional nature 
of the high temperature superconductors. In particular, 
we have neglected the difference between the transverse 
and tangential spatial directions in our treatment of the 
problem. The appropriate calculations would include this 
difference and lead to a non-symmetric shape. However, 
we neglect these complications at this stage. 

In order to make an analogy with the QCD case where 
the condensate on the core is described by a complex 
field, we are free to represent two degrees of freedom 
represented by a unit Neel vector rh, rh 2 = 1 defined 
by eq. (|19fl in terms of a single complex field $ as, 



in 



$ + $* 1-|$| 



i + |$| 2 '?:(i + |$| 2 )' 1 + 



(26) 



where $ = |$|e m (this is simply the projection of the 
unit sphere onto the complex plane) . At this point we are 
free to pick the direction of the Neel vector. For a back- 
ground classical field describing a vortex defined along 
the z direction, we will pick rh to lie in the xy-plane so 
that m z = 0. We should note that all calculations and re- 
sults which follow do not depend on the particular choice 
of m that we have made above. However, we do expect 
that this will turn out to be the lowest energy configu- 
ration when higher order derivative terms are included 
in the free energy (see below). We neglect fluctuations 
of the absolute value |$| for description of the classical 
background and consider variation of its phase a. In this 
case, we have |4>| = 1 for the classical background field as 
it follows from the transformation l|26(l , and rh simplifies 
to: 



1 



m 



($ + <FV(<r - $),o). 



(27) 
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with | $| = 1 to be fixed. The condensate |$| ^ can 
carry currents and charges along the string so we will 
represent it by the following ansatz which describes the 
dependence of these excitations on z, t: 

$ = \<S>\e la{z ' t] = \<f>\ e *( kz -^) . (28) 

With this redefinition of the fields, the kinetic term es- 
quires two additional terms due to the (z,t) dependence 
of the phase in the core: 



X 
2 



\d m\ 



|m 2 (r) [(d a(zM 



-v 2 \d z m\ 2 ] - 
v 2 (d z a(z,t)) 2 ] 



(29) 



where v 2 = p/x- The key point of the time dependent 
ansatz (|28ll is as follows. Naively, one could think that 
the time dependence in a classical solution brings an ad- 
ditional energy into the system which usually does not 
help to stabilize the configuration. However, as Witten 
noticed in if there is a conserved charge in the sys- 
tem, the configuration could be stable due to the con- 
servation of the corresponding charge. In a sense, the 
time-dependent configuration becomes the lowest energy 
state in the sector with a given non-zero charge. A similar 
time-dependent ansatz for a different problem was also 
discussed by Coleman in [35| where he introduced so- 
called Q-balls, macroscopically large stable objects with 
a time dependent wave function. We follow ref. 01 an d 
define a charge N which is topologically conserved: 



(30) 



where the path C is defined along the vortex loop and 
we assume that w, k are some constants along the loop. 
Since a can change by multiple of 27r in circling the vortex 
loop, N must be an integer. This is required in order for 
the condensate to to remain single valued. 

In addition to the topologically conserved winding 
number AT, there also exist the standard Noether charges 
and currents which can be trapped on the vortex core 
associated with the parameter u>, k included in the phase 
a above. In our case, the relevant symmetry, 50(3), 
implies conservation of three Noether charges: 

Qk = [d 3 rj° k = i X fd 3 r[(d m a ) (S k ) ah m b ], (31) 



while the corresponding three currents are: 

Jk = I d 3 rf k ,= -zp ( d 3 r[(d z m a ) (S k ) ab m b ], (32) 



where Sk are the three generators of SO (3). 

A vortex loop with nonzero Noether charges Q k and 
topological charge N trapped on the core is described 
by an ansatz of the following form which depends on 
the position f = {r,(f>,z) and time t as (using Eq.'s 
P )l .p)l. ipn) l. and ®): 



J4> 



(33) 



in = cry —g(r) (cos(/cz — ujt), sin(fcz — Lot), 0), (34) 

where as in the previous section /(r), g(r) are solutions to 
the classical equations of motion obeying the appropriate 
boundary conditions. 

For the solution given by (|34|l we have a nonzero 
Noether charge Q z which is trapped on the vortex core: 



Qz = xLtoT,, 



(35) 



where £ is defined as as the integral of \rh\ 2 over the 
vortex cross section: 



d 2 i 



(36) 



We should note that for a different choice of the direction 
of Neel vector to. the conserved charge which is nonzero 
would be different. The important point is that a nonzero 
charge will always be present independent of the of the 
Neel direction to. 

As we mentioned at the beginning of this section, these 
vortons are spinning and carry a non-zero angular mo- 
mentum. The vortons are stable against shrinking due to 
conservation of angular momentum. To calculate the an- 
gular momentum of a vorton with nonzero charges N, Q z 
trapped on the core, we use the standard formula for 
the angular momentum expressed in terms of the energy- 
momentum tensor: 
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d 3 r(T 0l Xj - T jXi), 



(37) 



which can be approximated for a large vorton in the plane 
as: 



M ~ 2ir X R 2 <J 2l -f I d 2 rg(r) 2 Lok, 
= 2nxR 2 u k S. 



(38) 



The angular momentum points in the direction normal to 
the surface formed by the vorton. We now see that the 
reason we would expect such configurations to be clas- 
sically stable is simple, it is just because these vortons 
are spinning and angular momentum is conserved. One 
can say that the vorton is stable because it is the low- 
est energy configuration in the given sector with nonzero 
conserved charges N,Q Z . Angular momentum M 138|l . 
which is essentially the product of two charges N and Q z 
is also nonzero when both charges N and Q z are non- 
zero. In the discussion above we neglected the higher 
order derivative terms. In particular, there will be some 
correlation between the charge Qi (|31|l and the momen- 
tum along the vorton Pi = T^i in the expression for the 
free energy (~ T^Qi). Such a correlation implies that 
the ansatz (|77|) will represent the lowest energy config- 
uration if the angular momentum M 138(1 points in the 
direction normal to the surface formed by the vorton. 

We will assume that we now have a vorton configura- 
tion with nonzero values of iV and Q z . In order to assign 
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specific numbers for these quantities one must look at the 
mechanism of formation. We will not address such com- 
plex issues in this paper and simply assume that there is 
some nonzero probability for a vorton to form. For recent 
work on the issue of vorton formation we refer the reader 
to [2f|. The free energy of a vorton can be obtained by 
substituting (J33J) and J3IJ| into (|2"U|) : 



T 



| ((uj 2 + v 2 n 2 )m 2 + v 2 (\7 rm ) 2 



+ v 2 \V^\ 2 ))~(\a\ + ~g)m 2 -\a\ 



—Mm 

2 v 



|V>| 4 + 2to 2 H 2 ) 



(39) 



where rn = tj^/\a\/b g(r). We can simplify this expres- 
sion further by using the fact that m is a solution to the 
equation of motion and represent the free energy in the 
following way: 



T = L hwrM ln(A/£) - ^£4 + X^^j , 
where we have defined the quantity E4 for brevity: 



(40) 



(41) 



The first term in (|4U|) is simply the energy from the dSC 
vortex with no condensate present in the core (to log- 
arithmic accuracy). Here A is the long distance cutoff 
which must be included to regulate the logarithmic di- 
vergence of the normal global string. The long distance 
cutoff is typically the distance between vortices, so in our 
case we will take A = L where L is the length the vortex 
loop. The second term is negative, reflecting the fact that 
it is energetically favorable to have an AF core. And the 
third term is the additional contribution to the energy 
due to nonzero Q z (N). 

There are various cases which must be considered, 
v s k > uj,v s k < w, and v s k = uj. Notice that the ef- 
fect of adding a having k,uj nonzero is the addition of a 
masslike term for m to the Lagrangian: 



As Davis and Shellard point out in [22] , if uj starts out less 
than v s k quenching occurs and forces uj increase faster 
than Q z /L. In the opposite case where uj > v s k anti- 
quenching occurs and therefore uj increases more slowly 
than Q z /L. The important conclusion that was drawn 
from this analysis is that uj/{v s k) — ^ 1 is an attractor 
[22| . As a loop shrinks uj/(v s k) approaches 1 and the 
quenching (or anti-quenching) slows and the eventually 
stops leaving a classically stable vorton behind. There- 
fore, for simplicity we will focus on the so called chiral 
case when uj — v s k which is the most stable configura- 
tion. Realistically, we know that the periodic structure 
of the material of the superconductor breaks the rota- 
tional symmetry, leading to non-conservation of angular 
momentum (it can be transferred to the material). How- 
ever, the topological charge N is still a conserved quan- 
tity. Therefore, for the chiral vortons, Q ~ N is also 
conserved due to the relations (|30|l and (|35|1 . For such 
configurations, stability is ensured. When uj ^ v s k, the 
transfer of angular momentum to the lattice is possible, 
eventually settling to the chiral case with uj — v s k. 

In the chiral case the size of the condensate £ is in- 
dependent of N,Q Z ,L and the free energy l|40|l can be 
written as: 



T = Li 



- (2n) 2 p— - 



(44) 



2 E4 ' 



where a s is the string tension of the bare dSC/AF vor- 
tex with Q z = N — and k is expressed in terms of 
the conserved winding number N according to eq. ij3U)l . 
Written in this form, it is immediately obvious that for a 
given nonzero value of N the free energy has a minimum 
at L = L Q : 



(45) 



We can give a crude estimate of the winding number 
density, no = jf- , of a stable vorton configuration: 



SC 



-{uj 2 



„2;„2 



k 2 )m 2 . 



(42) 



If v s k > uj then the effect of a nonzero v s k,u> is to add 
a positive mass term for \m\ to the Lagrangian. This 
counteracts the effects of the negative mass term in the 
original Lagrangian l|14|l . Since k ~ 1/L quenching oc- 
curs and the size of the condensate £ decreases as the 
vortex loop gets smaller. Conversely, if uj > v s k one has 
the opposite situation and anti-quenching occurs. As the 
vortex loop shrinks, the size of the condensate S gets 
larger as one would expect. The different cases have been 
examined using numerical calculations in [2(], • Recall 
that uj is given as 



(43) 



(46) 
(47) 



where 6 m is the width of the condensate. This gives us: 



n 



N 



1 



\g\ - y/*x(j* 2 - 



(48) 



which is approximately the inverse width of the conden- 
sate. As expected, the winding number density no does 
not depend on the large number N, but depends only on 
the internal structure of the vorton, i.e. on the width of 
AF condensate in dSC vortex core. The equation l|48|) 
tells us that as one goes around a vorton the direction of 
the Neel vector rh varies over a distance scale ~ <5 m , the 
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width of the condensate. As the doping is decreased and 
the AF-dSC phase boundary is approached from above 
the width of the condensate increases. For a given value 
of N(Q Z ) (determined at the time of formation) the size 
of a stable vorton increases as one approaches the AF- 
dSC phase boundary. 

The discussion above has shown that vortons are in- 
deed classically stable. This would imply that on the 
quantum mechanical level such quasiparticles are at least 
metastable. The issue of quantum stability of vortons 
was addressed in a recent paper [3jj. In this paper they 
calculated the lifetime of a vorton with Q — 0, N ^ 
(pure current case). The mechanism of decay is some 
quantum mechanical tunneling process where the con- 
densate instantaneously goes to zero on the core, allow- 
ing the winding number to decrease by one unit from N 
to TV - 1. 



IV. CONCLUSION AND FURTHER 
SPECULATIONS 

In this paper we have reviewed the dSC vortices which 
have an antiferromagnetic core within the 50(5) theory 
of high temperature superconductivity 0, EI LI ■ We 
have compared these dSC/AF vortices with similar vor- 
tices which arise in a completely different context, high 
density QCD OH OH, D3 . 

The main point that was presented in this paper is 
that loops of dSC/AF vortices called vortons can exist 
as classically stable objects in the presence of zero exter- 
nal magnetic field. The source of the stability of these 
vortons is conservation of angular momentum that coun- 
teracts the string tension, which prefers to minimize the 
length of the vortex loop. The fact that there is a conden- 
sate trapped on the vortex core is crucial for the stabil- 
ity of vortons. It is the condensate which allows nonzero 
charges to be trapped on the core, leading to the presence 
of nonzero angular momentum. It remains to be seen if 
such quasiparticles will be important for the physics of 
the high T c superconductivity. In what follows, we pro- 
vide arguments supporting the idea that the vortons can 
play a key role in AF-dSC phase transition. At this point 
we consider the vorton mechanism driving AF-dSC phase 
transition as a conjecture. 

The first argument goes as follows. As we have shown 
above, for a given value of N there exists classically 
stable vorton configurations with size L and fixed ra- 
tio N/Lq ~ 1/S m , where S m is the width of the con- 
densate that is trapped on the vortex core. As one de- 
creases the doping parameter and approaches the AF- 
dSC phase boundary, the width of the condensate 5 m in- 
creases. This is the direct consequence of the fact that the 
asymmetry parameter \g\ becomes smaller and smaller 
when the phase boundary is approached. From the rela- 
tion Lq/N ~ S m given above, this would imply that L , 
the length of a classically stable vorton, must increase. 
Decreasing the asymmetry parameter \g\ ~ 4%(/i 2 — [J? c ) 



further would result in a large vorton with large core 
size. The volume of the regions filled with the AF state 
behaves like LqS^ ~ |£f| -3 / 2 - When the phase transition 
line is approached, the regions with the AF state fill the 
entire space. At this point the AF-dSC phase transition 
occurs. Although our approximations are no longer valid 
at this point, because our description assumes that the 
vorton core size is much smaller than L and the interac- 
tion between strings can be neglected (plus many other 
assumptions not to be mentioned). These assumptions 
certainly fail in the vicinity of the phase transition. Nev- 
ertheless, the fact that the size of the AF regions inside 
of dSC phase increases rapidly when the AF-dSC phase 
transition is approached should be considered as a strong 
argument in favor of the vorton mechanism driving the 
AF-dSC phase transition. 

We would also like to make the observation that on 
the other side of the phase transition boundary, in the 
AF state, there are quasiparticles whose cores are in dSC 
phase |38j . Therefore, one can imagine a situation when 
one type of quasiparticles (dSC vortices with an AF core) 
becomes a different type of quasiparticles (AF skyrmions 
with a dSC core) when the doping parameter decreases 
and the phase transition line is crossed. 

The next natural question to ask is as follows: let us 
assume that vortons are indeed the relevant quasiparti- 
cles which drive AF-dSC phase transition at small tem- 
peratures. Can the same vortons be an essential part 
of the dynamics when the temperature (rather than the 
chemical potential fi) crosses the superconducting phase 
transition at T c ? If the answer is positive, we would have 
a nice unified picture for two different phase transitions 
on the (T, fi) plane. We believe the answer, indeed, could 
be positive (see arguments below). 

We start by reminding the reader that the pseudogap 
phase is characterized by the temperature T c < T < T* , 
when the Cooper pairs are already formed but the long- 
range phase coherence sets in only at the much lower 
temperature T c <C T* . It is believed that in this regime 
the phase order is destroyed by fluctuating vortices of 
the Cooper pair field ij) above T c j^. It is quite nat- 
ural to identify our vortons (loops of vortices) sliced 
by a two dimensional plane with vortex-antivortex pairs 
with distinct experimental signatures from ref. |39|. In 
this case, since underdoped cuprates are effectively two- 
dimensional, at finite temperature the loss of phase or- 
der may be expected to proceed via the Berezinsky- 
Kosterlitz-Thouless phase transition. In this case, the 
vortons discussed in the present paper, being sliced by 
the two-dimensional plane, become the vortex-antivortex 
pairs analyzed in ref. and could be responsible for 
the phase transition at T c . However, the picture of the 
phase transition here is quite different from what we pre- 
viously discussed regarding the AF-dSC separating line. 
In the present case, when T crosses T c the transition hap- 
pens because of the vortex-antivortex interaction which 
is proportional to pn^- hx(xi — X2) and not because the 
seeds of a new phase (the vorton cores) fill the entire 
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space. This is the typical two dimensional form due to 
the global nature of the vortices (local strings do not pos- 
sess this feature). The volume occupied by the vortex 
cores at this point is still much smaller than the volume 
of the system. It is well known that such a logarith- 
mic interaction is a key element for understanding the 
Berezinsky-Kosterlitz-Thouless phase transition. 

Encouraged by the argument given above, we extend 
our conjecture and assume that the same quasiparticles, 
vortons, are responsible not only for the AF-dSC phase 
transition but also for the phase transition separating 
the pseudogap and dSC phases at temperature T c when 
fi > fi c . The natural question to ask is: how does the 
critical temperature T c (fj,) depend on the chemical po- 
tential /i within this conjecture? To answer this ques- 
tion we recall that the critical temperature T c for the 
Berezinsky-Kosterlitz-Thouless phase transition is pro- 
portional to the strength of the logarithmic interaction 
mentioned above. In our case it is nothing but the string 
tension determined (mainly) by the vacuum expectation 
value of ip field, (|44|l . i.e. T c ~ a s . When /i = /i c 
the asymmetry parameter is zero and (ip 2 ) = \a\/b is 
determined by the unperturbed coefficients a, b. The 
simplest way to determine how the expectation value of 
the field i/j varies when [i increases is to introduce the 
asymmetry parameter in the "symmetric" manner |40j. 
Fasym = ~ g{rri 2 — \^P\ 2 ) such that negative g corresponds 
to the condensation of the \ij)\ 2 field, and positive g corre- 
sponds to the condensation of the \fh\ 2 field. This asym- 
metry parameter enters the string tension in dSC phase 
as follows, a s ~ (\a\ — g)/b, g < 0. From this expres- 
sion one can immediately deduce that a s (and therefore, 
T c ) will increase with /i when g being negative, becomes 
larger and larger in magnitude. It is interesting to note 
that the very same conclusion has been reached in the 
original paper Q , where the increasing of T c with chem- 
ical potential was explained as a result of mass increase 
of the 7r triplet (see [l| for details). 

Having presented our arguments supporting the con- 
jecture that vortons might be the relevant degrees of free- 
dom in the dSC phase at zero external magnetic field, 
we now conclude with a few remarks on how this pic- 
ture can be experimentally tested. First, the possible 
experimetal methods (such as fiSR and inelastic neutron 
scattering) for observing AF vortex cores were discussed 
in ref. 2] and we shall not repeat their analysis. However, 
we should mention that the fact that the AF cores do ap- 
pear in the vortices |EIE0I1;@ suggests that the 5*0(5) 
model of high T c superconductivity may be correct. Our 
original remark here is that AF vortex core size is on the 
order of dSC coherence length far away from the point 
g = 0, and becomes larger when the chemical potential 
approaches fj, c at zero external magnetic field. We expect 
that the average size (L) of vortons grows with temper- 
ature, and therefore, correlations between AF cores will 



grow with temperature as well. Similar behavior is ex- 
pected to occur when /i approaches fj, c at a fixed tem- 
perature. In this case the effect is expected to be even 
more pronounced because the volume occupied by the 
AF cores grows as |ff| -3 ^ 2 as discussed above, and there- 
fore the correlations as well as the magnitude of the local 
electron magnetic fields should scale accordingly. This 
picture suggests that the AF correlation length is propo- 
tional to L and could be very large, much larger than 
any other scale of the problem. Apparently, such large 
AF correlation lengths have already been observed in |(j 
and we would like to argue that this correlation is related 
to our vortons. One should remark here that such a large 
correlation length cannot be simply explained by the in- 
teraction between vortices (which have size £ ~ 20 A) 
because it would lead to a strong dependence on the Neel 
temperature as a function of the intervortex spacing con- 
trolled by the external magnetic field, while observations 
suggest that the Neel temperature is field independent 

All the effects previously mentioned require the ex- 
istence of the AF core in the vortex and they are not 
specifically sensitive to the existence of vortons, which is 
the subject of this work. The main feature of the vor- 
tons is that they can carry angular momentum (see eq. 
(|2Hl) and provide large AF correlation lengths (see dis- 
cussion above). Therefore, these excitations should be 
present if a dSC sample is rotated with nonzero angu- 
lar momentum. The situation is very similar to the 3 He 
and A He systems where superfluid vortices can be stud- 
ied by rotating liquid helium in a can. In many respects 
our vortons are similar to rotons, and presumably can be 
studied in a similar way using the technique developed 
for these systems. 

Finally, if vortons are indeed the relevant degrees of 
freedom in the high T c superconductors, it provides a 
unique opportunity to study cosmology and astrophysics 
by doing laboratory experiments in condensed matter 
physics. Over the last few years several experiments have 
been done to test ideas drawn from cosmology (see the 
review papers 0, for further details) . 
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